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Abstract. We present a briel pedagogical guided tour ol the most recent applications 
of nextensive statistical mechanics to well defined nonlinear dynamical systems, ranging 
from one-dimensional dissipative maps to many-body Hamiltonian systems. 

1 Introduction 

Thermodynamics and Boltzmann-Gibbs (BG) statistical mechanics constitute 
central formalisms of contemporary Physics. It is therefore of extreme impor- 
tance to clearly establish under what exact conditions they are expected to ap- 
ply, especially in what concerns the important and ubiquitous stationary state 
referred to as thermal equilibrium. This wisdom is more subtle than it looks at 
first sight. For instance, surprisingly enough, it is not yet exactly known the nec- 
essary and sufficient conditions for the celebrated BG distribution to correctly 
describe the thermostatistical state of a system (see, for instance, (!]]). In order 
to settle such conditions, it is necessary to turn onto the microscopic dynamics 
of the system [^|, would that be classical, quantum, relativistic, or any other one 
that appropriately describes the system under study. 

Some of the central ingredients of BG statistical mechanics are the well 
known logarithmic entropy (Sbg = — k Y^iLi Pi m general, and Sbg({Pi = 
1/W}) = khiW for equiprobability, where for simplicity we are referring here 
to the case of discrete distributions of probabilities), and its exponential weight 
Pi oc e~ l3Ei (fi = 1/kT, {Ei} being the energy spectrum) at equilibrium. This 
entropy has the remarkable property of extensivity for systems which are inde- 
pendent. More precisely, if we have a system composed by subsystems A and 
B such that p£ +B = pfpf (V(i,i)), then S BG {A + B) = S BG {A) + Sbg(B). 
This property is at the basis of the standard understanding of thermodynam- 
ics. Such a property is not necessarily universal because it cannot be justified 
on general arguments (see, for instance, (3)). Strong evidence exists nowadays 
that it should not be gencrically assumed, contrarily to a widely spread, text- 
book belief among not few physicists. Very specifically, the Boltzmann principle 
Sbg{{Pi = 1/W}) = klnW should not be used for any given specific system 
unless properly justified. It was Einstein Q the first who warned about this. 
However, many physicists (see, for instance, Q) have not yet done their way 
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through this important concept. To be more precise, let us focus on ./V-body clas- 
sical Hamiltonian e?-dimcnsional systems including two-body interactions which 
do not present any mathematical complexity (such as an infinitely attractive 
potential) at short distances. It is acquired by now that the usual BG thermal 
equilibrium correctly describes the t — ► oo stationary state for all systems whose 
interactions are short-ranged (e.g., if a two-body attractive potential decreases 
like l/r a with a/d > 1). In particular, in this case, it is irrelevant to consider first 
the t — > oo limit and then the N — > oo limit, or the other way around. The situa- 
tion is much more subtle if we have long-ranged interactions (e.g., < a/d < 1). 
In this case, the lirn/v^oo lim^oo still leads to the BG equilibrium, whereas the 
linif-joo limjv^oo generically does not. It happens, typically, that more than one 
basin of attraction exists in the space of the initial conditions. For some initial 
conditions the system directly goes to the BG statioanary state; for others it 
first goes to a non-BG state, remains there for a time which diverges with N, 
and only eventually goes to the BG state. In other words, depending on the size 
of the system and on its specific physical time scale, it might happen that the 
system remains in a quasi-stationary state different from the BG one for times 
longer than the age of the universe. In such cases, the only physically inter- 
esting situation is the non-BG one, which, perhaps for a vast class of systems, 
might be the one emerging within the nonextensive statistical mechanics we 
shall describe in the next Section. The present paper is a tutorial review of these 
concepts as they naturally emerge for three simple but important classes of sys- 
tems, namely, (i) one-dimensional dissipative maps (e.g., the logistic map[^J), 
(ii) two-dimensional conservative maps (e.g., the standard mapQ), and (iii) 
classical many-body long-range-interacting Hamiltonians (e.g., the HMF model 



[p|Jg| J10| , p. 1|]12| and one of its generalizations, the so called a — XY model 13,14]). 
We shall not address here the many applications of nonextensive statistical me- 
chanics available in the literature. The interested reader may refer to Jl5| for 
various reviews. Just as an illustration, let us mention some of those applica- 
tions: turbulence jl^,^7), electron-positron annihilatio n |l8| , diffusion of Hydra 
viridissima fl9(| , diffusion of quarks in gluon plasma | 2(fl, Levy and correlated 
anomalous diffusions |2lj] , linguistics p2f , economics p5 , 24 1 , fluxes of cosmic rays 
[ Bp , solar neutrinos |26|] , high energy particle collisions [|27f , self-organized criti- 
cality among others. In one way or another, these phenomena seem to share 
long-range correlations in space/time, either long-range microscopic interactions, 
or long-range microscopic memory (nonmarkovian processes), or (multi) fractal 
boundary conditions, or, generically speaking, some mechanism which creates a 
scale-invariant hierarchical structure of some sort. Such possibility appears to 
emerge in classical many-body Hamiltonians and other nonlinear dynamical sys- 
tems, cverytime the Lyapunov spectrum approaches zero, hence chaos becomes 
impossible in the sense that the sensitivity to the initial conditions diverges less 
than exponentially in time, though a weak mixing can remain. 
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2 Mathematical formalism 

The first and main step in order to go from BG to nonextcnsive thcrmostatistics 
is to propose the use of an entropic form which generalizes that of BG, as follows 
& 

Sq = k ^±E±- (£ Pi = l;g6R). (1) 

9 i=l 

(for simplicity, and without loss of generality, we shall adopt k = 1 from now 
on). This nonnegative form can be conveniently rewritten as 

S q = (In, -> , (2) 

Pi 

where (...) = ^■_ 1 (...)pj, and the q-logarithm function is defined as 

- 1 

hiq = — (lnj x = ln.T; x > 0) . (3) 

The inverse function is the q- exponential one, given by 

e° q = [l + (l-q) x ]ih (e^ = e x ), (4) 

solution of dy/dx = y q with y(0) = 1. S q satisfies several remarkable properties, 
such as concavity for all {pi} (Vq > 0), stability with regard to {pi} (Vg > 0) f30| , 
and, for A and B independent subsystems, the following pseudo-extensivity: 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) , (5) 

from which the denomination nonextensive comes. For completeness, it is conve- 
nient to mention at this point that similar entropies are the Renyi one, defined 
as = [hiX^l^iPf]/[l — 9b an d the normalized one = (intro- 
duced independently by Abe and Rajagopal |n| and by Landsberg and Vedral 
| f?2| . Renyi entropy is very useful in the context of multifractal analysis (of non- 
linear chaotic systems, for example). However, S q and are not concave for all 
positive values of q, nor are they stable |3(J. These facts constitute severe draw- 
backs for basing a thcrmodynamical formalism on them. It is however convenient 
to have in mind that, since they are related to each other through monotonic 
functions S* = ]n[l + (l-q)S q }/[l-q] and = S q /[l + (l-q)S g ], they become 
extremized by the same probability distributions (assuming the supplementary 
constraints for optimization to be the same). 

3 Applications to low-dimensional maps 

Low-dimensional systems are a field in which the nonextensive statistical me- 
chanics has found one of its first and most proficuous applications. What makes 
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low-dimensional systems extremely attractive is the simplicity of their dynam- 
ics and the possibility to have direct access to the phase space, sometimes even 
to plot a visual image of the entire phase space. In fact, the iterative rules 
of low-dimensional maps can be easily implemented on a computer and have 
provided useful methods to perform "numerical experiments" to check for the 
validity of the nonextensive formalism. As mentioned in the previous section, 
the BG statistics is expected to fail when the sensitivity to initial conditions is 
not exponential. Therefore, different authors have looked and found applications 
of the nonextensive formalism to cases in which the maps have a nonstandard 
sensitivity to initial conditions, as indicated by the vanishing of the Lyapunov 
exponents. Different techniques have been produced in order to calculate the 'en- 
tropic parameter' q, and it turns out that this parameter characterizes 'classes 
of universality' of systems, at least for some of their statistical behaviors. For 
pedagogical reasons here we will discuss only some of the results obtained for 
the logistic map, as representative of a dissipative system, and for the standard 
map, as an example of a conservative one. 

3.1 Dissipative Maps: the logistic map 

The logistic map is one of the simplest one-dimensional dissipative system one 
can imagine. It has been intensively studied and it has led to important de- 
velopments in chaos theory [^|. Since it has all fundamental characteristics of 
non-conservative systems, the logistic map is often taken as a paradigmatic ex- 
ample; in other words it can be considered as the Ising model of non-linear 
dynamical systems. The logistic map can be described by the following iterative 
rule]] 

x t+ i = fa(x t ) = l-ax 2 t - 1 < x t < 1; < a < 2; (6) 

and it shows different regimes according to the value of the control parameter a 
[HIHl • F° r small values of a the map is regular since it has non-positive Lyapunov 
exponents (see Fig.Q). In particular, for a < | the iterates of the logistic map 
converge to a single fixed point x* = f a (x*), whose value depends on a. At a = | 
we have A = and we observe the first pitchfork bifurcation: for a value slightly 
larger than | the iterates converge to a cycle-2. The number of points of the 
attractor f a (x) keeps doubling at distinct, increasing values of the parameter a, 
until a critical value a c = 1.401155198... is reached, where the attractor becomes 
infinite. This is known in the literature as the period-doubling route to chaos. In 
fact, beyond a c the behavior of the logistic map is chaotic (positive Lyapunov 
exponent A) for most of the values of a. At the chaos threshold a = a c , A is once 
again zero, and this point is the well known edge of chaos which has been the 
subject of intensive debates about the definition of what a complex systems is 

1 An equivalent, perhaps more popular version of the logistic map is x t +i = rx t (l — x t ), 
with < x t < 1 and 1 < r < 4. However, we notice that form (ph is easily generalized 
to other universality classes of order z > 1: x t +i = 1 — a|a;t| z . In this way, a z-logistic 
map is in fact representative of the universality class of unimodal maps of order z; 
i.e., of one-dimensional maps with a single maximum of order z. 
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Fig. 1. (a) Attractor of the logistic map as a function of a. The edge of chaos is at 
the critical value a c — 1.401155198... (b) Lyapunov exponent A as a function of a. 



P'4p'I • I' i s exactly at the edge of chaos and at all other points where A is zero 
(as at period doublings, but also for instance at inverse tangent bifurcations |^|), 
that a generalization of the standard BG entropy will be fruitful, as we show 
below. 
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Sensitivity to initial conditions Let us start by explaining better the mean- 
ing of the crucial concept of sensitivity to initial conditions. For all the cases in 
which the Lyapunov exponent A is positive (negative) we expect on the average 
an exponential increase (decrease) of any small initial distance £ (t) = J*'-x' j *° 
hold. We typically have 

£(i)=exp(At), (7) 

where x t and x£ are the positions at time t of two initially close trajectories that 
may be considered, for the sake of generality, in a e?-dimcnsional phase space. 
In the case of the logistic map d = 1. This case is referred in the literature as 
strongly sensitive (insensitive) to initial conditions. What is not usually reported 
in textbooks is what happens when A, as defined by equation (^), is equal to 
zero. In Ref. |36) was proposed that in such situations the system still exhibits 
sensitivity to initial conditions, though in a form that is the g-gcneralization of 
Eq. ©: 

^(t)=eiq ?q (X q t) = [l + (l-q)X q i\Th ( ? gl), (8) 

This equation recovers Eq. (Q) for q = 1, and displays a quite rich spectrum of 
possiblities according to the sign of \ and q. Let us concentrate on the case 
X q > 0, q < 1. In this case the system presents a weak sensitivity to initial 
conditions, in the sense that two initially close trajectories diverge in time as a 
power-law instead than exponentially. A weakly sensitive system is well described 
once q and X q are given. In |36| q was calculated numerically for the logistic map 
at the edge of chaos, observing that the upper bound of £(i) lies on a line that is 
precisely a power- law (see Fig. ||); furthermore, it was conjectured that a value 
q = 0.2445.... could be deduced from the Feigenbaum's constant a. These results 
were recently confirmed in J37J, where it was proved that, at the edge of chaos 
the upper bound of £(t) has exactly the form of Eq. (||), with q = 0.2445... 
and Xq = lna/ln2 = 1.3236.... In fact, as a consequence of the Feigenbaum's 
renormalization group recursion relation, two close initial conditions starting in 
the neighborhood of the origin x = 0, produce, each t = 2™ — 1 (n = 0, 1, 2, ...), 
a dominant power-law separation of their iterates that has the form (|^). After 
these particular steps, the iterates reproduce self-similar sequences, that are also 
power-laws. The effect of starting the initial points in other regions of the phase 
space instead that in the vicinity of x ~ is just a shifting in time of Fig. |^. 

The renormalization group approach has proven powerful also at other points 
where the Lyapunov exponent is zero. In ]38[ ] the same q-exponential form (||) 
has been exactly obtained for the sensitivity to initial conditions at pitchfork and 
inverse tangent bifurcations, now with other values for q and X q . It is important 
to notice that the previous scheme of application of Eq. (||) generalizes to the 
other universality classes of unimodal maps of ordcr^] z > 1. 

Entropy production The analysis of the entropy production provides an al- 
ternative approach to the same problem and a second way to characterize the 



See previous footnote. 
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Fig. 2. Sensitivity to initial conditions for the logistic map at a — a c . The dots represent 
for two initial data started at xo — and x' ~ 1CP 8 . The solid line is the function 
in Eq. (|), with q = 0.2445... and A, = In a/ In 2 = 1.3236.... 



g-index at the edge of chaos. Instead of considering two trajectories initially close, 
we take a distribution of initial conditions localized in a tiny region of the phase 
space, we let this distribution evolve in time according to Eq. (|fij), and we study 
the increase of the entropy |3j|§| . In the case of the logistic map, we can partition 
the phase space interval — 1 < x < 1 into W equal cells, and consider an initial 
distribution of N points placed at random inside one of the cells; the initial cell 
is also chosen randomly in the partition. The normalized number of points that 
occupy a cell defines a probability distribution: pi = rii/N = 1)j an d 

as the system evolves, from the probabilities Pi{t) we can calculate the entropy 
production S q (t), via Eq. (|j). In the chaotic regime the BG entropy (q = 1) 
exhibits, before a saturation due to the finiteness of W, a linear increase in time, 
and in this stage the entropy production rate (that is expected to be equal to 
the Kolmogorov-Sinai entropy pCfl ) coincides with the positive Lyapunov expo- 
nent This is shown in Fig7~|| for the case a = 2, where S q (t) is reported 
for three different values of q. As t evolves, S q (t) increases (in all cases bounded 
by W \lq X i or hiW^ when q = 1), but only the curve for q = 1 shows a clear 
linear behavior with a slope equal to the Lyapunov exponent A = In 2 . For 
q < 1 the curve is convex, while for q > 1 the curve is concave. The slope in the 
linear stage does not depend on the dimension of the cells of the partition and on 
the size of the initial distribution. Therefore, in the chaotic regime the standard 
B-G entropy is the only g-entropy that displays a finite entropy production rate 




Fig. 3. Time evolution of S q for the logistic map with a = 2. We consider three different 
values of q. The results are averages over 500 runs with N = 10 6 and W = 10 5 . 



limj—too w^oo S q (t)/t. We wish now to repeat the same analysis at the edge of 
chaos, where A is zero. Due to the power-law sensitivity to initial conditions and 
to the fractality of phase space (the attractor now is a multifractal) we expect 
that anomalous behavior may be observed. In Fig. |4| we consider a = a c and we 
plot S q (t) for four different values of q; the curves are obtained with a large W 
and with an average over many different runs since the complexity of the phase 
space originates many fluctuations. Consistently with the value of q extracted 
for the sensibility to initial conditions, the result is that the growth of S q (t) is 
found to be linear when q = 0.2445..., while for q < q c (q > q c ) the curve is 
convex (concave). This behavior is similar to the one of Fig. ||, with a major 
difference: the linear growth is not at q = 1, but at q = 0.2445... Q. 



Multifractal analysis There is also a third, geometrical method that allows 
us to calculate the q of the logistic map. This method is based only on the de- 
scription of the multifractal attractor existing at a = a c , and gives exactly the 
same value we have already obtained. The multifractal attractor can be char- 
acterized by using the multifractal function /(a) |33|. This function is defined 
in the interval [a m i n , a max ], and its maximum equals the fractal or Hausdorff 
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Fig. 4. Time evolution of for the logistic map with a = a c . We consider four different 
values of q. The case q = 1 is in the inset. W = N = 2.5 10 6 and results are averages 
over 15115 runs. 



dimension dr. The value of q can be calculated from 



1 1 



1 — q a 



rain Q-raax 



(9) 



where a rn i n {a m ax) are the two values for which f(a) is defined. In particular 
for the logistic map at the edge of chaos a m in = 0.380..., a max = 0.755... and 
we have again q = 0.2445.... 



3.2 Conservative Maps: the standard map 

Conservative systems present quite different statistical properties than dissi- 
pative ones. One of these differences is that the transition from chaoticity to 
regularity happens in the phase space without an 'edge of chaos', but in a rather 
complicated way as described by the Kolmogorov-Arnold-Moser (KAM) theo- 
rem see Fig. ||e-g). Conservative maps may be derived, for example, from 
Hamiltonian systems. In this case the original system is characterized by an even 
number of dimensions d — 2n, where n is the number of degrees of freedom. If 
the Hamiltonian H is time-independent the dimension can be reduced by one. 
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Fig. 5. (a)-(d) Phase portrait of the standard map for typical values of a. N = 20 x 20 
orbits (black dots) were started with a uniform distribution in the unit square and 
traced for < t < 200. (e)-(g) Islands-around-islands. (e): N = 100 x 100 initial data 
were started inside a small square of side ~ 10 -2 and traced for < t < 5000. (f) is 
a magnification of the island inside the rectangle in (e). (g) is a magnification of the 
island inside the rectangle in (f). 



Moreover, as in statistical mechanics we are interested mostly in recurrent tra- 
jectories (i.e., in those trajectories that come back again and again, indefinitely, 
to any part of the phase space they have once visited), we can take a Poincare 
section of the phase space cutting transversally the constant-energy hypersurface 
and considering the successive intersections of each orbit with this transversal 
surface. In this way we reduce the dimension of the phase space to d,M = 2n — 2 
and we obtain a great numerical simplification, as we can now analyze a discrete- 
time iteration map instead of a continuous-time system of differential equations. 
The map thus obtained may be shown to be simplectic (see, e.g. [^2[); this im- 
plies that the cIm Lyapunov exponents are coupled in pairs, where each member 
of the pair has the opposite of the other. 

The minimum dimension that allows for the application of this scheme is 
oIm = 2 (n = 4) , and a paradigmatic example of this kind, that plays for con- 
servative systems a similar role the logistic map plays for dissipative ones, is the 
standard map, that may be characterized by the equations^: 



xt+i — yt + — sin(27ra; t ) + x t (mod 1), 

Z7T 



Vt+i = Vt + TT" sm ( 27ra; t) (mod 1), 

Z7T 



(10) 



3 The map is also known as the kicked rotator map. 
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where ogi (( = 0,1,..,), The system is integrable when a = 0, while, roughly 
speaking, chaoticity increases when \a\ increases (see Fig. ||a-d). 

As for the logistic map, we expect that anomalous effects may be produced as 
the complexity of the phase space increases (i.e., when we pass from the chaotic 
to the regular regime reducing \a\), and as this happens, we would like to see 
how the nonextensive formalism generalizes the BG one, at least for some of the 
statistical variables we have previously described. 



Sensitivity to initial conditions In the case of the standard map, the sen- 
sitivity to initial condition depends dramatically on the position in the phase 
space. Chaotic regions produces exponential separation of initially close trajec- 
tories, while regular regions exhibit linear separation of initially close orbits. 
At the border between these regions, appear fractal-like structures of islands- 
around-islands (see, e.g., 43 1, see Fig. ||e-g). We notice however that in statisti- 



cal mechanics we are mainly interested in extracting global averages behaviors, 
so that we can try to calculate an average £ (t) by sampling uniformily the whole 
phase space. When this experiment is performed for rather small values of |a|, as 
we show in Fig. |6|, a crossover in time between two different regimes occur J44| . 
Fig. ^a displays in fact a linear increase of In (£(£)) only after a certain iteration 
step. Correspondingly, Fig. ^b exhibits a linear increase of ln 9 (£(i)) precisely 
for those steps where ln(£(i)) is concave (numerically it was found q ~ 0.3). 
Moreover, the crossover time increases when \a\ decreases. This means that two 
different behaviors dominates the average sensitivity to initial conditions for dif- 
ferent time. Initially, because of the predominance of fractal-like structures in 
the phase space, the sensitivity to initial conditions is power-law; then, after a 
certain time, the sensitivity to initial conditions becomes exponential due to the 
rapidity of the exponential growth if compared with the power-law. If we wish to 
model such a behavior with a single function, we can first observe that exponen- 
tial and power-law sensitivity to initial conditions may be viewed as solution of 
the differential equations £ — A£ and £ = A g £ 9 respectively. A differential equa- 
tion describing a crossover in time is then: £ = A£+(A 9 — A)£ 9 (q < 1, A, A g > 0). 
The solution [Eta 



£(*) = 



Xq (l-g)At 
A + A 



(11) 



presents, in the case < A <C A g , three asymptotic behaviors, namely (i) linear: 

i i i 

£ ~ 1+Xqt, for < t < t crossl = {1 _ q)Xq ; (ii) power-law: £ ~ [(1 - q)X q ] i 1 -*, 

for t crossl «t« tcrossi = (x- q )\ 'i ( ni ) exponential: £ ~ (x) 1 ^^*: for 1 > 

tcross2 • 



Entropy production If the connection between sensitivity to initial condition 
and entropy production that we have observed for the logistic map is correct, 
we should observe a crossover between different statistical regimes even when 
we analyze the entropy production of the standard map. This has been done 
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Fig. 6. Average sensitivity to initial conditions of the standard map for small values 
of a. |xo — Xq ~ 10~ 9 and results are averages over 5000 runs. An initial separation of 
order In (a) ln(£(t)) displays a linear increase only after a certain time. Correspondingly, 
in (b) lno. 3(£(i)) exhibits a linear increase for the initial steps. 

in Efj], with a numerical experiment analogous to the one previously described 
for the logistic map and that is difficult to perform because of the complexity 
of the phase space and the related effects of rapid saturation of the entropy for 
W finite. In Figs. and || we display how, for small values of \a\ and before 
saturation begins, the entropy production of S\(t) becomes linear only after a 
stage where S q (t) is linear, with, once again, a value q ~ 0.3. Notice that the 
crossover time tends to infinity when \a\ goes to zero, so that in this limit the 
stage characterized by the power-law sensitivity to initial conditions and by this 
nonstandard entropy production lasts for infinite time. 

4 Applications to Many-particle systems 

4.1 Hamiltonian nonextensive systems: the HMF model 

In the last years there has been an intense research activity in investigating 
microscopic dynamical features in connection with macroscopic features. In par- 
ticular several authors have focused their attention in characterizing from a 
dynamical point of view macroscopic phase transitions |4l| . [Is[[l9H50]]5l| . It has 
been found, for example, that in correspondence of a phase transition there is 
a peak in the chaoticity of the microscopic dynamics. In general chaoticity is 
of fundamental importance for any system to reach the standard equilibration. 
Anomalies should be expected when this does not occur, since ergodicity break- 
ing might happen and the system can remain trapped in some region of phase 
space. In this respect a model whose features have been particularly interesting 
is the Hamiltonian Mean Field model, from here on referred as HMF model, 
and its generalizations where the range of the interaction can be changed. The 
model was introduced in ref || and extensively studied in refs ]^,p|, |lO| , pl][r^JT-| . 
In this volume one chapter is dedicated to its equilibrium properties both from 
a dynamical and from a thermodynamical point of view J12[. In the present 
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Fig. 7. For the standard map we show Si(t) for different values of a < 2 (5000 to 7000 
runs were averaged). Full circles correspond to JV = W = 1000 x 1000 for a = 2, 1.5, 1.1 
(N = W = 5000 x 5000 for a = 0.65); empty circles correspond to N = W = 448 X 448 
for a = 2, 1.5, 1.1 (N = W = 2236 x 2236 for a = 0.65). Inset: Slopes of Si(t) (ifi(flO) 
and of S . 3 (t) =20 (if .3 ^ 0.71 + 0.30 a 3 6 ) in their linear regimes (see also Fig. |). 
The lines are guides to the eye. 

section we consider only the anomalies found in the out-of-cquilibrium regime 
before equilibration to the Boltzmann-Gibbs regime and how these behavior can 
be related to nonextensive statistics . 

The HMF model describes a system of N planar classical spins interacting 
through an infinite-range potential^. The Hamiltonian can be written as: 



where 9i is the ith angle and pi the conjugate variable representing the angular 
momentum or the rotational velocity since unit mass is assumed. The interaction 
is the same as in the ferromagnetic XY model [||, though the summation is 
extended to all couples of spins and not restricted to first neighbors. Following 
the Kac procedure |T4| , the coupling constant in the potential is divided by N. 
This makes H only formally extensive (V ~ N when N — > oo ) |2^ . [l5| . p^|Jl^1 , since 
the energy remains non-additive, i.e. the system cannot be trivially divided in 
two independent sub-systems. The canonical analytical solution of the model 
predicts a second-order phase transition. At low energy density the system is in 
a ferromagnetic phase characterized by a magnetization M ~ l ; where M is the 



N „.2 1 N 
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Fig. 8. For the standard map we show So.3(i) for: a — 1.5 (a); a = 1.1 (b); a = 0.9 (c); 
a = 0.7 (d); a = 0.65 (e) (5000 to 7000 runs were averaged). Full circles correspond 
to N = W = 1000 x 1000 for (a), (b), (c), N = W = 5000 x 5000 for (d), (e); empty 
circles correspond to N = W = 448 x 448 for (a), (b), (c), N = W = 2236 x 2236 for 
(d), (e). Inset: t croas (a) defined as the intersection of the linear part (before it starts 
bending) and a standard extrapolation of the bended part of the curves 5*0.3(4). 
Notice that t C ross(a) diverges for a — > 0. The lines are guides to the eye. Our results 
suggest limi^oo lim a ^o limvp-»cx> limjv^oo s °-^ ~ 0.71 for q = q* ~ 0.3, whereas this 
limit vanishes (diverges) for q > q* (q < q*). 



modulus of 

1 N 

M=-^m, , (13) 

i=i 

with m, = [cos(8i), sin(9i)]. Increasing the energy beyond a critical value, the 
spins become homogeneously oriented on the unit circle and M ~ 0. The caloric 
curve, i.e. the dependence of the energy density U = E/N on the temperature T, 
is given by U = ^ + ^ (1 — M 2 ) P^l anc ^ snown m Fig|] as a full curve. The crit- 
ical point is at energy density U c = 0.75 corresponding to a critical temperature 
T c = 0.5 ||. The dynamical behavior of HMF can be investigated in the micro- 
canonical ensemble by starting the system with water bag initial conditions, i.e. 
9i = for all i (M = 1) and velocities uniformly distributed, and integrating 
numerically the equations of motion || . As shown in Fig.|^, microcanonical sim- 
ulations are in general in good agreement with the canonical ensemble, except 
for a region below U c , where it has also been found a dynamics characterized by 
Levy walks, anomalous diffusion and a negative specific heat|ll[]. Ensemble 
inequivalence and negative specific heat have also been found in self-gravitating 



systems |52|, nuclei and atomic clusters 53 54j.p5|.p6| . In our case such anoma- 



lies emerge as dynamical transient features, although in a generalized version of 



the HMF model ml 12 it seems that they remain also in the BG equilibrium 
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u 

Fig. 9. We compare the equilibrium caloric curve of the HMF model (full curve) with 
the numerical simulation for a system of size N=500. The numerical integration is 
followed up to a short time, when the system is in the metastable quasi-stationary 
state (QSS). 



regime. In order to understand better this disagreement we focus on a particular 
energy value, namely U = 0.69 below the critical energy, and we study the time 
evolution of temperature, magnetization, and velocity distributions. 

In Fig.[l(] we report the temporal evolution of 2 < K > /N for U = 0.69 
and N = 500. This quantity, at equilibrium, is expected to coincide with the 
temperature (< • > denotes time averages). The figure shows that when the 
system is started with out-of-cquilibrium initial conditions (water bag in our 
case) rapidly reaches a metastable quasi-stationary state (QSS) which does not 
coincide with the canonical prediction. After a short transient time, not reported 
here, 2 < K > /N shows a plateau corresponding to a N-dependent temperature 
Tqss(N) (and Mqss ~ 0, M — ► if N — > oo) lower than the canonical temper- 
ature. This metastable state needs a long time to relax to the canonical equi- 
librium state with temperature T can = 0.476 and magnetization M can = 0.307. 
In the same figure we plot a typical single-particle motion in the two insets. In 
panel (b) the motion is characterized by anomalous diffusion and Levy walks 
with constant velocity [j^. This motion is well described by a model introduced 
in ref. |p7| . These features change at equilibration, see panel (c) where the mo- 
tion is a standard random walk and diffusion becomes normal. In ref. [fiof it 
has been found that the lifetime of the metastable QSS and the crossover time 
from anomalous to normal diffusion do coincide within the numerical accuracy 
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Fig. 10. For the HMF model at U=0.69 and N=500, we plot the temperature, cal- 
culated by means of the average kinetic energy, as a function of time. Two plateaux 
emerge from the calculation. The first one is the metastable quasi-stationary state 
(QSS), while the second one corresponds to the Boltzmann-Gibbs (BG) equilibrium 
regime. In the inset (b) we show the angle vs time of a typical rotor in the metastable 
regime: the motion shows Levy walks and anomalous diffusion . In the inset (c) we 
show the single-particle motion in the BG equilibrium regime: the motion is a random 
walk. 



adopted. On the other hand, the duration of the plateau increases with the 
size of the system. In fact the lifetime of QSS has a linear dependence on N as 
shown in Figjll] where simulations for different out-of-equilibrium initial condi- 
tions are plotted. These numerical results indicate that the two limits t — > oo 
and N — > oo do not commute: if the thermodynamic limit is performed be- 
fore the infinite time limit, the system does not relax to the BG equilibrium 
and the anomalies discussed above remain forever . Indeed similar features have 
also been found for spin glasses [ p8|j59| ] and granular matter Q,^0| and could 
probably be interpreted within a general theoretical scenario. When N increases 
Tqss{N) tends to a definite value which is = 0.38. The latter is obtained 
analytically as the metastable prolongation, at energies just below U c = 0.75, of 
the high-energy branch charaterized by M = 0. In rcf. |^l| we have also found 
that [Tqss(N) — Too] oc TV" 1 / 3 . This scaling, exploiting the expression for the 
caloric curves, implies also a scaling of the magnetization, which goes to zero like 
Mqss TV" 1 / 6 . A further numerical check of this law was obtained in ref . |6^| . 

It is important at this point to discuss the behavior of the Lyapunov expo- 
nents as a function of N and of time in order to investigate if also this measure 
shows some anomaly. In fact this is just the case: also the largest Lyapunov 
exponent A has two temporal regimes and in particular two plateaux. We plot 
in Fig.|l2| the time evolution of A for the case U = 0.69 and N = 500. The mi- 
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Fig. 11. For the HMF model we plot the lifetime of the metastable quasi-stationary 
state, tqss, for U=0.69, as a function of N. Different out-of-equilibrium initial condi- 
tions are indicated by different symbols. The figure shows that the duration of these 
states diverges with N. 



mcrical simulation averaged over 350 events shows that in the metastable QSS 
regime the largest Lyapunov exponent is smaller than at the BG equilibrium. In 
the inset we report the scaling of A in the QSS regime as a function of the size: 
increasing N, Xqss tends to zero as 

N -x/9 H Thc scaling law for thc Lya _ 

punov exponent can be related to the scaling of the magnetization as discussed 
in ref. j62|. Thc fact that A — > for N — > oo implies that mixing is negligible 
as N increases and one expects anomalies in the relaxation process in the sense 
advanced by Krylov ]63] | and in a very similar way to the case of the maps 
discussed in the previous section. 

In order to check relaxation in HMF not only the temperature, the second 
moment, but also the higher moments of the velocity distribution should coincide 
with those of the BG equilibrium one. In Fig.|l3| we focus our attention on the ve- 
locity probability distribution functions (pdfs). The figure shows that the initial 
uniform velocity pdfs, quickly acquire and mantain for the complete duration of 
thc metastable state a non- Gaussian shape . In particular the velocity pdf of thc 
QSS is wider than a Gaussian (drawn as dashed curve) for small velocities, but 
shows a faster decrease for p > 1.2. The enhancement for velocities around p ~ 1 
is consistent with the anomalous diffusion and the Levy walks (with average ve- 
locity p ~ 1) observed in the QSS regime flfj|| . The following rapid decrease for 
p > 1.2 is due to conservation of total energy. The stability of the QSS velocity 
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Fig. 12. We plot the largest Lyapunov exponent, A, vs time for the HMF model at 
U — 0.69 and TV = 500. Two plateaux are found also in this case. In the inset we show 
that Xqss goes to zero as TV _1//9 . This behavior can be explained with the anomalous 
fluctuations of the magnetization. See text for further details. 



pdf can be explained by the fact that, since for TV — > oo, Mqss — > also the 
force on the spins tends to zero with N, being Fi = —M x sin9i + M y cos6i. If 
N is finite, one has only a small random force, which makes the system even- 
tually evolve into the usual Maxwell-Boltzmann distribution after some time. 
When this happens, Levy walks disappear and anomalous diffusion leaves place 
to Brownian diffusion |Tc|], see insets in Fig.[l(]. Nonextensive thermodynamics is 
able to reproduce the non-Gaussian pdf reported in Fig.^. In fact the nonex- 
tensive formalism provides, for the canonical ensemble, a q-dependent power-law 
distribution in the variables pi , Qi. This distribution has to be integrated over 
all 9i and all but one pi in order to obtain the one- momentum pdf, P q (p), to 
be compared with the numerical one, Pnumip), obtained by considering, within 
the present molecular dynamical frame, increasingly large N-sized subsystems of 
an increasingly large M-system. Within the M » TV >> 1 numerical limit, we 
expect to go from the microcanonical ensemble to the canonical one (the cut-off 
is then expected to gradually disappear as indeed occurs in the usual short-range 
Hamiltonians), thus justifying the comparison between P q {p) and P n um(p)- The 
enormous complexity of this procedure made us to turn instead onto a naive, 
but tractable, comparison, namely that of our present numerical results with the 
following one-frce-particlc pdf HHI P (p) °c [1 - (wX 1 ~ q)p 2 } 1/{1 ~ q) , which 
recovers the Maxwell-Boltzmann distribution for q = 1. The same type of for- 
mula has been recently used to describe successfully turbulent Couette- Taylor 
flow and non-Gaussian pdfs related to anomalous diffusion of Hydra cells 
in cellular aggregates. In the case of HMF the best fit is obtained by a curve 
with entropic index q = 7, and temperature T = 0.38. The agreement between 
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Fig. 13. We plot the velocity probability distribution functions of the HMF model 
for U=0.69 and different N values (open symbols). The Gaussian equilibrium shape is 
also plotted (dashed curve) together with the fit with the nonextensive pdf formula 
obtained for q — 7 and T = 0.38 (full curve). 



numerical results and theoretical curve improves with the size of the system. A 
finite-size scaling confirming the validity of the fit was reported in rcf. J6l| . Since 
in our case we get q > 3, the theoretical curve does not have a finite integral 
and thus it needs to be truncated with a sharp cut-off in order to make the total 
probability equal to one. It is however clear that, the fitting value q = 7 is only 
an effective nonextensive entropic index. In ref. [pi] we verified, through the cal- 
culation of the fractal dimension Di , that a dynamical correlation emerges in the 
/i-space before the final arrival to a quasi-uniform distribution. During intermedi- 
ate times some filamentary structures appear, a similar feature has recently been 
found also in self-gravitating systems |30), which might be closely related to the 
plateaux observed in Fig.|l(| We have found that these correlations do not sensi- 
bly depend on N, thus likely the possible connection does not concern the entire 
^i-spacc, but perhaps only the small sticky regions between the "chaotic sea" and 
the quasi-orbits |65|| . Recently another application of the generalized nonexten- 
sive formalism has been observed for this model in the relaxation regime. It has 
been found that the time relaxation of macroscopic quantities, like the temper- 
ature for example, obey power-laws and more specifically q-cxponential curves. 
To illustrate this, in Fig.|lJ, we plot the time relaxation of the temperature to 
its final value for U = 0.69 and for two sizes N=500,1000. In particular we show 




Fig. 14. For th HMF model, we plot the difference between the temperature T and 
that one at equilibrium T eq , for U=0.69 and N=500,1000, vs time (open symbols). The 
points are the result of an ensemble average over 1000 runs. The figure shows that 
the relaxation process is a power law in time and can be very well reproduced with 
a decreasing g-exponential. We report the g-exponential fits for comparison (full and 
dashed curves). In the inset we plot the q- logarithm of the same data. The entropic 
index q is related to relaxation and should not necessarily be equal to that one of 
velocity pdfs. 



y(t) = \T — T e q\, with T eq = 0.476 as a function of the time t. The system relaxes 
to equilibrium according to the q-exponcntial formulan 



y(t) = 2/o 



i-(i-s) 



t l V(l-9) 



r 



(14) 



where q is the entropic index, r a characteristic time and yo a saturation value 
[ p9[ . For N = 500 the same data shown in Fig.JT^ are used. For this particular 
size, one finds a value of the entropic index equal to q — 1.71 and a characteristic 
relaxation time r = 18000. A slightly greater value of q equal to 1.81 is found for 
N=1000. In the inset we report the q-logarithm of same data (6(|. This entropic 
index q is related to the relaxation and there is no reason to be exactly the same 
as the one obtained for the velocity pdfs. It is interesting to notice that similar 
relaxation process characterized by a g-exponential has been recently discovered 
also in echo experiment on the Internet |j67f . Two final considerations before 
concluding this section. The first one is that in rcf.Bq] it has been found that 



4 The minus in front of the term (1 — q) at variance with formula (Q) is due to the fact 
we are considering a decreasing g-exponential and q > 1. 
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the HMF model shows aging in correspondence of the metastable QSS. Thus 
it is a very intriguing challenge left for future investigation, the study of the 
possible connections between spin glasses, hamiltonian mean field models and 
nonextensive statistics. Second, the anomalies we have found are not a peculiar 
characteristic of the HMF model. In fact the same behavior has been found also 
in a generalization of the HMF model, the a — XY model jD|[l4j, in which the 
spins are placed on a d-dimensional lattice, and the interaction energy between 
two generic spins i and j (of formula (|l2|)) is modulated by an extra factor 
Here Tij is the distance between spin i and spin j If ajd < 1 nonextensivity 
induces anomalies similar to those described for the HMF model. Conversely, if 
a/d > 1, i.e. when the potential decreases rapidly enough with r, the a — XY 
model behaves in the standard way and does not exibit any of the anomalies 
discussed in this section. Similar features induced by nonextensivity have been 
found also for long-range Lennard- Jones- like potentials |39| . 

5 Conclusions 

We have finally ended our guided tour in the realm of nonextensive systems. 
Our choice was to focus only on a few simple examples, chosen in order to 
represent the most instructive classes of systems. On one hand we have consid- 
ered low- dimensional systems, in particular low-dimensional maps as the logistic 
map (representative of a dissipative system), and the standard map (represen- 
tative of conservative systems). On the other hand we have considered high- 
dimensional systems discussing the anomalies observed in the HMF model, a 
classical many-body long-range-interacting Hamiltonian. Though representing 
very different kind of systems, all the examples reported share some common 
properties: long-range correlations in space/time, either long-range microscopic 
interactions or long-range microscopic memory (nonmarkovian processes), or 
(multi)fractal boundary conditions, or, generically speaking, some mechanism 
which creates a scale- invariant hierarchical structure of some sort. Such proper- 
ties appear in classical many-body Hamiltonians and low-dimensional nonlinear 
dynamical systems everytime that the Lyapunov exponents approach zero, i.e. 
when chaos (strong mixing) is impossible and the sensitivity to the initial con- 
ditions is not exponential in time. In most cases of weak sensitivity to initial 
conditions, a generalized nonextensive formalism may replace standard thermo- 
statistics. 
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